
Institute for Space Studies 





l EVOLUTION OF O STARS. III. HELIUM BURNING 


■;> 



Richard Stothers 


Gpo 

cfsti 


price $ 
PRICE(S) $ 



Hard C °P V (HC) 
M 'crofiche ( MF ) 

ft 6 53 July gg 



M 66-12865 

H fACCESSlON NUMBER) 

QA 


(PAGES) j 




(THRU) 


-X 


'5x0 


(CATEGORY) 


Goddard space Flight Center 
National aeronautics and Space Administration 


k 



1 



HELIUM BURNING 


Richard Stothers 


Institute for Space Studies 
Goddard Space Flight Center, NASA 
New Yorfc, N.Y. 



2 


ABSTRACT 

izfir 

The evolution of a star of 30 M is considered during 

0 

the phases of helium ignition (I), depletion (He), and 
exhaustion (X) . 

After the brief I-phase (1.5 x 10^ years) during which 
the onset of helium burning halts the gravitational contrac- 
tion in the core, helium is depleted in a large convective 
region growing outward in mass fraction until it nearly 
overtakes the hydrogen burning shell. Helium depletion 
causes a radius shrinkage until Y c = 0.3, whereupon the 
radius re-expands. The final C /O ratio in the core is 

uncertain, but will be zero if the critical reduced alpha 
16 

width in O ' is greater than 0.3. The luminosity of the 

hydrogen shell source is low and monotonically decreases, even 

though the shell finally reaches the hydrogen-rich semi- 

convective zone (aX = 0.15) when Y =0.01. The star 

c 

covers the spectral range from Bl to B8, while semi- 
convection maintains the total luminosity at a nearly 
constant value. The lifetime of the He-phase is 5.2 x 10 5 


years 
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4 

During the X-phase (1.3 x 10 years), the carbon-alpha 

reaction is able to sustain the luminosity in the shrinking 

-4 

convective core until Y c = 10 . Thereupon the gravita- 

tional contraction in the core dominates the nuclear 

3 

sources, and rapidly (10 years) brings the star into the 
region of late-type supergiants. The weak shell burning, 
prolonged by hydrogen enrichment due to full convective 
mixing at the inner boundary of the semiconvective zone, 
will then be extinguished by the envelope expansion. Since 
the new convective zone (along with the semiconvective zone) 
eventually retreats toward the surface, leaving behind a 
radiative zone of frozen composition, two major fixed 
composition discontinuities occur in the stars hydrogen- 
helium (q = 0.44) and helium-oxygen (q = 0.37). 

The evolutionary features of very massive stars are 
discussed generally, and a time scale computed for the very 
late phases with and without neutrino emission. In either 
case, the total lifetime to iron core formation will be 
slightly less than 6 million years. 
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I. INTRODUCTION 

The evolution of massive stars during the advanced 
phase of helium burning is not well known. Among the stars 
with M > 10 M^, this phase has been studied only in the case 

of 15 M (Iben 1964) and 15.6 M (Hayashi, Jugaku, and 

o 0 

Nishida 1959 f 1960; Hayashi and Cameron 1962; Hayashi, Hoshi, 
and Sugimoto 1962) . The present paper considers in detail 

the evolution of a very massive star of 30 M from helium 

o 

ignition through exhaustion of central helium, and sketches 
briefly the more advanced phases up to iron core formation. 

It is a continuation of calculations on the previous phases 
(Papers I and II, Stothere 1963b, 1964). 

II. ASSUMPTIONS 

We adopt the same general assumptions of Papers I and 
II, including the complete mixing of material in convective 
regions and the absence of overshooting at the boundaries 
of convective and semiconvective regions. The distribution 
of energy sources is more complicated, however, as will be 
discussed in the next section. In conformity with the pre- 
vious calculations, the opacity is assumed to be due only to 
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electron scattering, and the equation of state is represented 
by the sum of the perfect-gas and radiation pressure. The 
core is found to remain non-degenerate. 

The age- zero chemical composition, retained by the 
outer radiative envelope, is again taken to be 

Xe- 0.70 , X- 0 . 2.7 , Ze- 0.03, * C M0 = (1) 


III. GENERAL STRUCTURE 

The overall structure of a very massive star is 
different in each of the evolutionary phases considered in 
this paper. We isolate three phases, as follows. 

1. Helium ignition . — As the gravitational contraction 
of the core following central hydrogen exhaustion proceeds 
to raise the core temperature, helium burning begins near 
the center. The mass of the tiny convective region stops 
shrinking and increases again. Very quickly helium burning 
halts the gravitational contraction entirely, but the 
hydrogen shell source is still required to maintain the total 
luminosity. 
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2. Helium depletion , — During this phase helium is 
depleted in a large convective core growing out in mass 
fraction. Since the time scale is long compared with the 
previous phase, hydrogen depletion in the shell is also 
significant, and the mass fraction of the shell peak increases. 
Toward the end of central helium burning the shell reaches 

the semiconvective zone and burns at its inner boundary. 

3. Helium exhaustion . — As the helium content of the 
convective core drops below 1 per cent, gravitational con- 
traction again begins to play an important role. The 
contraction is further accelerated because of incipient 
photoneutrino emission. Having attained a maximum mass, the 
convective core now begins to retreat and leaves behind a 
region composed chiefly of oxygen. The hydrogen shell 
source continues to burn very weakly. As the total luminosity 
increases, the innermost part of the semiconvective zone 
becomes fully convective . 

We adopt the following system of notations to designate 
the stellar zones, as closely in accordance with the notations 
of Papers I and II as possible: (I) the homogeneous outer 



7 


radiative envelope retaining the age-zero chemical composi- 
tion, (Ila) the semiconvective intermediate zone, (lib) the 
convective intermediate zone, (IIla,b,c) the radiative inter- 
mediate zones, (IV) the convective core. Zones Ilia, Illb, 
and IIIc are separated by the hydrogen-helium and helium- 
oxygen discontinuities, respectively. The following sub- 
scripts are also adopted: 1 for the Interface I-II, 2 for 

the inner boundary of Zone Ila, 3 for the outer boundary of 
Zone III, and 4 for the Interface III-IV. The hydrogen- 
helium discontinuity (hydrogen burning shell) will be 
designated by a subscript s. Thus for a semiconvective shell, 
q g = q 2 = q^, whereas for a convective shell, q g = q 3 » The 
helium-oxygen discontinuity will be designated by a sub- 
script a* Thus for the outward-growing convective core, 

q = q Jt Quantities referring to the immediately preceding 

or 4 

evolutionary model will be marked with an t. 

The successive evolutionary phases will be denoted: 

I (helium ignition) , He (helium depletion) , and X (helium 
exhaustion) , followed by a number indicating the particular 
model in the phase. Figure 1 represents the stellar 
structure in the phases considered in this paper. 
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Values of solar luminosity and radius were taken from 
Allen (1963), whereas in Papers I and II Chandrasekhar's 
(1939) values were used. 


IV. ENERGY RATES 
a) Nuc lear Energy 

Throughout all phases considered in this paper, the 
temperature of the hydrogen burning shell does not differ 
significantly from 5 x 10 ° K, so that we may use the 

following interpolation formula for the rate of the CNO-cycles 

£ h - XXu.0 fT*, (2) 

O 

with v = 14 and log g = -99.0. We assume the full equili- 

H 

brium abundance of oxygen. 

12 16 

Below the hydrogen shell, C and 0 will certainly 

have reached their equilibrium abundances by the end of the 

hydrogen burning phase. Hence X /X « 0.96 or X « 0.015. 

N CNO N 

By the onset of helium burning most of the will have 

, 14 4 ip 

been destroyed by the reaction N + He 0 , except for 
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14 

a small amount of N mixed down as the convective core moves 
outward. Consequently, the initial abundances in the core 
for the main phase of helium burning are assumed to be 

Y^O.^7, Xcn s O, Z -- 0*03 , (3) 

and at first the nuclear energy generation of the core will 

be due to the triple-alpha reaction. For simplicity, we 

also neglect the contribution from nitrogen burning in the 

helium ignition phase. 

12 

When C is formed by the triple-alpha reaction, the 

12 4 16 

C reacts with He to form 0 . Subsequent alpha particle 

captures by 0 and 0 produce small amounts of neon and 

magnesium. However, from the calculated models, it may be 

estimated (Reeves 1964) that the final abundances will be 

20 24 

less than about 5 per cent for Ne and Mg . Consequently, 
we are justified in considering the energy generation as 
due chiefly to the helium and carbon reactions. Reeves 
(1965) gives the energy generation rates as 
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He a 


X 10 


-i «V T « 




3.f 5- X 10° | Y 3 p* Tg 3 


(4) 
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e c = x io 17 |X c YpTj 51 (l- o.ooh-Tj to^T 9 ' 2/1 ) 2 x io" 8 , (5) 


where 


$= 30.0ST T 8 ( / 3 f o.Ofe T/ /? (\ - O.OG7 T s ‘ /3 ) y (6) 

for the reactions 3 He^ — > and + He^ — * 0^, 

respectively. Here T n = T X 10 -n . The weak electron 
screening correction applies at the densities of our models, 
and is given by (for both reactions) 

4- 1+ J.H- f* T 4 ' 5/> . < 7 > 

For simplicity we assume that f = 1. Afterwards, with the 
help of the calculated models, our assumption may be justi- 
fied by showing that f is close to unity. 

2 

The reduced alpha width, 0 , of the 7.12 Mev level in 

a 

16 2 
0 is still very uncertain. We have used 0^ = 0.1 with a 

possible error of a factor 10. Fowler and Hoyle (1964) give 

arguments for a rather large value of 0.73. We shall later 

o 

investigate the effects of altering the value of 0 , which 

a 

enters the rate as a straight multiplicative factor. 
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b) Gravitational Energy 

The rate of energy release due to gravitational con- 
traction in the core is given by (Paper II) 

[|c Xrv ) (8) 

where y = ( 1— 0 ) / 8 and the time derivative applies to a 
given mass fraction. 


_c) Neutrino Emission 

During helium exhaustion, neutrino losses begin to 
become important in dissipating energy from the star. The 
only process operative in our low-density models will be 
photoneutrino emission, y + ® -* e — + v + "y. The rate is given 

by (e.g. Reeves 1965) 

= S-XIO 7 T,® . 


(9) 
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V. BASIC EQUATIONS 

The basic equations for the stellar structure are 
given by equations (5) - (15) of Paper I. However, during 
the present phases L(r) is not constant outside the convec- 
tive core. Therefore in equations (11) we replace the 
envelope parameter C by C' = Cw, where w = L(r)/L. The 
formulas for the various sources of luminosity are given 
below. 


a) Shell Luminosity 


Let the energy generation rate be represented by 

equation (2) , and consider the shell source thin, with a 

peak at q and with X = 0 below q . Analogously to the 
s s 

procedure of Hayashi and Cameron (1962), we write 


^ % + a* % Z > 


( 10 ) 



/ v ¥ 
f sTs 


( 11 ) 



p ~V~ <£ JLocj ^ . 

u TIcff'J 


s 


( 12 ) 
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Now we may integrate p to obtain 

H 


L « ’ (-R + pt *- pf) 


(13) 


In the radiative intermediate hydrogen zone, the composition 

gradient is frozen, so that the coefficients of equation (10) 

remain constant. In the semiconvective zone the composition 

gradient is continually changing, but is much less steep 

than in the radiative zone; therefore we set x = X = a 

s o 

with little loss of accuracy. In the convective intermediate 

zone, X = X exactly. The exponent y is evaluated at q : 

3 s 


Y s X * 

3L[ 
u l 

J-J3 

. ru-»-| (3 

V 

f -JL 

, 1 “I 

y s u i 


'W q J 

V- ¥- 

■ t 


Y IT 


a J 


(radiative) (14) 


(semiconvective) (15) 


(convective) (16) 


where \ = d log u/d log q = (5/4) a (a q + 2a q 2 ) . Beyond 

1 2 

q g , we assume no energy sources, so that here w = 1. 
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where v will be a function of the central temperature, 
except in the case of the photoneutrino loss. Then, generally, 


L * J 

O 

°> *■ / (Yr tV'R*" 3 " l *,*/ Iff/ l t'» ) 


In particular, we have the following cases: 


L = 

L>He ) 

L c, 

L„, 


0) - 


1, 

o, 

( 19 ) 

K.' 


£ cl*c^, 

t 0 y . 
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c) Gravitational Luminosity 

Energy release due to gravitational contraction is 
assumed to occur below the hydrogen-helium discontinuity 
where the hydrogen burning shell is located. This assumption 
is later justified by the detailed calculations. From the 
results of Paper II the rate of gravitational energy release 
is approximately independent of mass fraction. Consequently 
we have 


^3 * e 3 ^ ~ % + J > 


(I-phase) (20) 


^ %9 


(X-phase) (21) 


In the case of the I-phase , it is found that the steep 
temperature dependence of the helium burning rate allows 
the nuclear source to sustain the convective core against 
its gravitational contraction. Thus we assume that gravi- 
tational energy release occurs only in the dehydrogenized, 
non-helium-burning Zone Illb. For the X-phase, is 
chosen to be q or q , as discussed in Section XI. 

S 



VI. INTEGRATION METHOD 


With trial envelope parameters (to be specified below) , 

Zones I - IV are integrated from the surface inward to the 

point in Zone IV where the conditions U < 3 and dU/dq < 0 

are no longer fulfilled. Revision of the trial parameters 

is made until these conditions are met at some prechosen 

mass fraction close to the center. Then a starting value 

for e is extrapolated, so that outward and inward inte- 
c 

grations may be performed. Fitting is made in U and V at a 
specified 8^ just inside Zone IV. 


VII. NUCLEOSYNTHESIS 


As hydrogen is consumed in a shell, the mass fraction 
of the shell moves outward at a rate given by 


Ax 


L-h 

ME h > 


( 22 ) 


where E 

H 


6.0 x 10 


18 


erg/gm and 
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In the helium burning core, the time rates of change 
of He 4 , C 12 , and 0 16 are given by 


££ * - ins. _ _l 

iz E He 3 E c 

(24) 

m 6r 

* 

(25) 

>Xo >*• tr 
* 3 Ec ' 

(26) 

Consequently, we may write for the core 


- A fo-il , L He 1 

Atr " KE Hc , 3 MEc 1 

(27) 

AjXc^ L„_ Lr 

At ME He M£ c ’ 

(28) 

X 0 = 1-Ze-TC-Xc . 

(29) 


Combination of equations (27) and (28) yields a quartic 


equation for Y. 
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VIII. UNSTABLE INTERMEDIATE ZONES 
a) Semi convection 

In Paper I it was shown how the following criterion 
for convective instability. 




> 




*- 3(5 


( 30 ) 


must be used in a chemically inhomogeneous region. The 
equality sign in this criterion is used to determine the 
distribution of chemical composition that is actually set 
up by mixing in an unstable zone, so that the zone attains 
a state of convective neutrality. This mixing phenomenon 
is called semiconvection. 

Since the upper boundary of the semiconvective zone 

(q^) merges smoothly into the outer radiative zone, the 

material here is nearly homogeneous; consequently, any 

convective instability is uninhibited at q . At the lower 

boundary (q^ = q ) the steepening gradient of mean molecular 

weight eventually chokes off the convective motions. Thus, 

in the early models during the hydrogen-burning phase, the 

zone below q 0 is actually radiative even though (n+1) >(n+l) 

2 ad rad 
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These considerations remove the apparent contradiction found 
by Savedoff and van Dyck (1959), who neglected the term in 
<3-tn|a/d^nP in equation (30) . 

The effect of semiconvection on the stellar evolution 
was shown to be small during the phases considered in Papers 
I and II. In fact, the effect is completely negligible 
throughout the hydrogen burning phase, since the following 
conditions are both satisfied at the same mass fraction 



- vj/v, * on 

-C X II H * I, 1 , *UI H ■ 02) 


Hence for computational purposes, it was sufficient to 
consider the intermediate zone of inhomogeneous composition 
as wholly radiative. Eventually, however, the growing semi- 
convective zone makes itself felt by slightly changing the 
values of the physical variables in the inner radiative zone 
(but not the chemical composition which remains frozen) from 
those values that would obtain if the semiconvective zone 
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were assumed radiative. This effect shows up when condi- 
tions (31) and (32) are not simultaneously or individually 
satisfied. By the onset of helium burning, condition (32) 
is satisfied at a mass fraction that is several per cent 
smaller than Computationally, we now integrate the 

semiconvective zone directly from q^, applying the jump in 
chemical composition (eq. [31]) when condition (32) is met. 
The jump in |j to the frozen value of |i in the radiative zone 
automatically produces the proportionate jump in U and V, 
This was not the case before, since and depended on 
the integration of the assumed radiative zone. 

b) Full Convection 

When the inner boundary of the semiconvective zone 
reaches the dehydrogenized helium zone, condition (32) is 
replaced by the condition that hydrogen be conserved: 

Aj*Xdq = 0. As the luminosity continues to increase, it 
eventually becomes impossible to satisfy this condition 
beyond a certain q^, denoting the maximum inward extent of 
the semiconvection. A new mechanism must be found to 
maintain a downward transport of hydrogen. Apparently the 
only selfconsistent scheme involves the growth of a fully 
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convective region between the seraiconvective zone and 
the radiative helium zone. In this scheme the semi- 
convective zone retreats with sufficient helium enrichment 
to maintain convective neutrality, while the growing 
convective region steadily drains hydrogen and helium from 
its inner boundary. The net effect is to increase the 
hydrogen content of the homogeneous convective region, since 
the semiconvective zone is always relatively richer in 
hydrogen. 

With modest increases in the luminosity, the inner 
boundary of the convective region moves very slightly inward 
of q m . As a good approximation, we assume composition dis- 
continuities at the inner (q^) and outer (q^) boundaries of 
the region. At these boundaries it is necessary that 

(n+1) = (n+1) , in order to preserve continuity of the 

rad ad 

radiative flux. This condition of equality, in fact, 

determines the amount of composition jump at q^ t since 

(n+1) < (n+1) for q < q < q . Finally, we require 

rad ad 2 1 

that hydrogen be just conserved when (n+1)^^ = ^ n+ ^^ a< a 

at q . 

3 

By the end of helium exhaustion, the inner boundary of 
the convective region begins to move outward. It is found 


I 



that a new radiative zone is left behind. Hence, both the 


semiconvective and fully convective regions start to 

approach the stellar surface. Within the new radiative 

zone left behind, the composition gradient is determined by 

t = l (q/q ) \ analogous to the gradient left behind by 
s s 

the shrinking convective core during hydrogen burning 

(Paper I) . However, in the present case, X can be directly 

determined by use of l and l . Thus for either an inward- 

3 m 

or outward-moving convective region, q^ is the only unknown, 
which must be solved for by trial and error until the 
condition of hydrogen conservation is met. 


IX. HELIUM IGNITION PHASE 

In this phase we adopt a run of chemical composition 
throughout the star as shown in Table 1. The hydrogen 
abundance in the radiative intermediate Zone Ilia is given 
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This is the distribution found in the last model of the 
gravitational contraction phase following hydrogen exhaus- 
tion (Paper II) . However, we now assume that X = 0 below 
an infinitely thin shell source at q^. The effect of an 
abrupt change of X at q through k and n is not large if X 
is small. Since the time scale of the helium ignition phase 
is short, we further assume that hydrogen depletion in the 

shell source is negligible; thus q will remain constant. 

s 

From Paper II, the peak of the shell source is at q = 0.328. 

s 

It is found that the models during the brief I -phase of 

unchanging chemical composition are not very sensitive to 

the particular value of q g ; therefore we have chosen the 

shell peak of Model G6, at which X = 0.018. 

s 

Integration of models for this phase is performed as 
follows. First, the evolutionary step is taken by specifying 
C. Since gravitational contraction is assumed to occur only 
in Zone Illb, the non-dimensional luminosity in this zone is 
given by 




( 34 ) 
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where w* = w /(q -q ). Trial integrations are then performed 
g g s 4 

with a fixed value of w*, for various values of w until a 

9 s 

fitted model is obtained. The fitted model yields q and 

4 

hence w . Substitution of the non-dimensional variables of 
9 

Paper I (eqs. [81 and [10]) into equation (13) gives a 
relation between w g and R, which then determines the physical 
variables. Since no significant depletion of the helium 
occurs, nuclear burning in the core is due only to the triple- 
alpha reaction. Hence the nuclear core luminosity may be 
calculated from equation (18) , yielding an implied value for 
the gravitational luminosity, w = 1 - w -w . A new 

g s He 

I 

value of w is obtained from w and q , and the whole proce- 
g g 4 

I 

dure is iterated until w = w . 

g g 

The time step, £t, is computed from equation (18) by 

using equation (20) to obtain g and by evaluating the 

y 

physical variables at a point midway between q and the mean 

s 

value of the present and previous q . For the I-phase only, 

4 

the semiconvective zone was treated in the manner described 
in Paper I . 

Table 2 gives the results for six models representing 
the helium ignition phase. The first model corresponds to 
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a stage only slightly more advanced than the last stage in 
Paper II. After a brief time, the steep temperature de- 
pendence of the helium burning rate (v « 20) allows helium 
burning to sustain the entire dehydrogenized zone by only a 

very slight increase in T and p . The convective core 

c c 

grows extensively in mass fraction as a consequence. Although 

the computed time steps are rather rough, it is estimated 

that hydrogen depletion in the shell, corresponds to an 

equivalent change in q of only 0.01. As we have noted 

s 

above, the structure will not be sensitive to q ; therefore 

s 

neither will be at. The total time elapsed between the 
onset of helium burning at the center and cessation of the 

4 

core gravitational contraction is about 1.5 x 10 years. 


X. HELIUM DEPLETION PHASE 

During the main phase of helium burning, the only 
energy sources are nuclear: helium and carbon burning in a 
convective core and hydrogen burning in a thin shell. The 
distribution of chemical composition is given by Table 1 
and equation (33) . Since the time scale of the He-phase is 
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long, nuclear burning in the shell depletes a significant 
amount of hydrogen and causes q to increase (Fig. 2<a) . 

S 

In the convective cone, helium burns into carbon and oxygen. 

However, the resulting composition discontinuity at the 

boundary of the outward-moving core will not violate the 

condition of continuity of radiative flux at non-nuclear 

burning interfaces (Paper I), since the electron scattering 

opacity is constant in a dehydrogenized region. 

In the numerical integrations, the evolutionary step 

is taken by specifying q . With fixed values of C and 

s 

(Y, X ) at the center, trial integrations are run with 

v 

various values of w until a fitted model is obtained. As 

s 

in the I-phase, equation (13) determines R and thus all the 
physical variables. The core luminosity is then obtained 
from equation (18) and gives an implied value for the shell 

I I 

luminosity, w = 1 - w - w . Comparison of w and w 

S H6 C S S 

yields an improved guess for C. The procedure is iterated 

I 

until w = w for some value of C. Then the time step, At, 
s s 

is computed from equation (22) , using the amount of hydrogen 

depletion, AjXdq, between q g and q g . The choice of q g for 

the initial model is not critical, because X is small over 

s 
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a large range in mass fraction near q = 0.328. Improved 

s 

values of (Y, X^) are obtained by solving equations (27) 
and (28) . Finally, the whole scheme may be iterated to 
obtain the definitive values of C and (Y, X ) . 

For the He-phase, the semiconvective zone was treated 

explicitly, as described in Section VIII. Its effect may 

be assessed by comparing the initial model for this phase 

(Table 3) with the final model for the previous phase 

(Table 2), where the semiconvective zone was assumed 

radiative. Both models represent the same stage. The 

effect is chiefly a slight lowering of the luminosity and 

a relative increase in the central condensation. The lower 

luminosity is due to the indirect effect of the hydrogen 

discontinuity at q This discontinuity steepens the 

—0 4 

radiative temperature gradient, dT/dr ~ hh ~ (1+X) * , 

beyond what it would be in the assumed radiative case for 
the same luminosity. Consequently, the shell temperature 
becomes too high, causing ^ T g v to run away. The shell 
temperature may be reduced only by decreasing dT/dr, hence 
by shrinking the hydrogen discontinuity (and therefore the 
extent of semiconvection) . Consequently, the total 
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luminosity must be lower in the semiconvective case. On 
the other hand, the relative increase in central condensa- 
tion, U ^ ~ yx \ is caused chiefly by the hydrogen enrich- 
ment of the inner envelope caused by semiconvection. 

When the shell source consumes its way to the semi- 
convective zone, AX is only 0.15 at the interface (Fig. 2si) . 
Thus the sole effect of the encounter is a slight prolonga- 
tion of shell burning, by which remains temporarily 
constant. The reason for the convective instability just 
outside the shell lies essentially in the steepening 
temperature gradient there. As pointed out in Paper II, 
(n+1) approaches (n+1) .as p decreases. However, in 

3.Q IT 3.0 S 

the present case, the decrease in p is due to the increase 

s 

in h ~ 1+X as the shell moves outward, since the total 

luminosity remains nearly constant during helium burning. 

This may be seen by observing that (1 -p)/hL ~ (n+l) ra( g « 

(n+1) . constant, 

ad 

For a wide range of stellar masses, the radius begins 

to expand when Y =0.3 (cf. 4 M and 15.6 M in Hayashi, 

c © © 

Hoshi, and Sugimoto 1962). Since the nuclear burning rate 

3 

goes as Y for helium burning, as opposed to X for hydrogen 
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burning, the value Y c = 0.3 corresponds to that value 

X c = 0.03 at which the radius changes direction during 

hydrogen burning. When Y c is as low as 1 per cent, the 

star still has a radius characteristic of an early-type 

supergiant. The reason is essentially that the nuclear 

core luminosity is enabled to continue increasing by carbon- 

alpha burning at high temperature (see Sec. XII). 

In order to maintain the high luminosity, £ c remains 

low. We note that 0 c is lower at the end of helium burning 

than at the end of hydrogen burning. This produces the 

more extensive convective core in the present case. By 

Model HelO, the core size has reached its maximum extent, 

only 0.07 in mass fraction below the hydrogen shell source 

(Fig. 2b). The low value of g also keeps the ratio p /T^ 

c c c 

small. Since nuclear reactions not only sustain the core 

but expand it slightly during the early stages of helium 

burning, p tends to remain very low. Consequently, T 
C O 

is high in order to maintain the appropriate nuclear 
reaction rate. It is because of these relatively higher 
temperatures in the cores of massive stars that the carbon- 
alpha reaction (which occurs at a higher temperature than 



30 


the triple-alpha process) is able to get under way before 

all the helium is burned into carbon. 

The products of the carbon-alpha reaction depend 

16 

critically on an uncertain reduced alpha width in O , as 

discussed in Section IV. With our adopted value of 
2 

0^ = 0.1, the core is still composed of one-third unburned 

carbon (by mass) at the end of helium burning. Calculations 

based on the model parameters of Table 3 indicate that the 

composition of the core will be all oxygen and heavier 
2 

elements if 9 a exceeds 0.3. 

XI. HELIUM EXHAUSTION PHASE 

In this phase, gravitational contraction in the core 
gradually replaces helium burning as the source of energy, 
when the abundance of helium drops below 0.01 by weight. 

The core composition does not change significantly, and in 
view of the uncertainty of the final products, we have 
simply adopted (a = 2 below q = 0.371, the maximum extent 
of the convective core (Model HelO) . The helium disconti- 
nuity at q = 0.371 (Fig. 2b) represents the location of 
the helium shell, which is not ignited during the stages 
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considered here. Hydrogen depletion in the outer shell 

is negligible on account of the short time scale involved; 

therefore q is determined by the convective instability, 
s 

as discussed in Section VIII. Table 1 illustrates the run 

of chemical composition throughout the star. It does not 

change appreciably during this phase. 

Since the gravitational contraction rate is not 

uniform between the composition discontinuities (q < q < q ) , 

a s 

we have constructed two sequences of models, with q^ = q^ 
and q^ = q^, as extreme cases bracketing the true models. 

Model construction proceeds as in the helium ignition phase 
(Sec. IX), except that we assume here gravitational contrac- 
tion occurring uniformly throughout the whole region below 

q and a central nuclear burning due only to the carbon-alpha 
9 

reaction (minus photoneutrino emission) . Hence we substitute 
W g = W g^g ^ nto e( 3 uation (34) and compute w^ =1 + w v - w g - w c « 
The time step, At, may in the present case be evaluated at 
the center of the star, and our guessed values for (Y, X ) 
in the convective core improved, as for the helium depletion 
phase. A sufficiently accurate expression for the time 
scale may be approximated from equation (8) s 
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3 _ lTT c 

^ ~ x i* c H' 


( 35 ) 


The neglect of A(i c in the early stages will not be serious 

since & is small. Furthermore, o remains nearly constant 
g M c 

so that we may also neglect Ay^. 

The results for the X-phase are displayed in Table 4. 
The assumption of q g = q g or q^ = q^ leads to significantly 
different temperatures and fractional luminosities at the 
tabulated stages. Consequently, only rough values of the 
radius and effective temperature are given. 

During the onset of the X-phase the region just out- 
side the hydrogen discontinuity at the shell becomes fully 
convective, mixing down additional hydrogen and prolonging 
shell burning (Fig. 2a). The boundary at q = q remains 

S O 

remarkably constant for moderate changes in the luminosity, 
although there is indication of a slight decrease. Even- 
tually, however, when the increasing central condensation 
necessitates an even higher luminosity, the convective 
intermediate region moves out as a whole toward the surface, 
continuing to penetrate the semiconvective zone but leaving 
behind a new radiative zone of frozen composition (Fig. 2a). 
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The inner boundary of this new Zone Ilia is thus fixed at 

q = 0.439 with X = 0.528. 
s s 

When g is stabilized as the gravitational luminosity 
c 

increases, the convective core shrinks in mass fraction. 

The final products of nuclear burning are deposited in a 
radiative zone left behind by the retreating convective core. 
As we have seen, this nuclear ash will be predominantly 
carbon and/or oxygen. The exhaustion of nuclear fuel takes 
place on a gravitational time scale, similar to that com- 

4 

puted for the I-phase, about 1 x 10 years. Of this, only 

3 

10 years are required to cross the Hertzsprung Gap. It 
is uncertain whether helium will be completely exhausted 
from the inner core before or after this crossing; however, 
nuclear burning in the core is found to become negligible 
when T c = 3.7 x 10 °K and = 2800 gm/cm . At this point 
photoneutrino emission will contribute 3 per cent of the 
luminosity. 

The internal evolution of the star is displayed in 
Figure 1, from the onset of helium burning through its 


exhaustion at the center. 
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XII. EVOLUTIONARY FEATURES 

In this section a general exposition of the basic 
features in the evolution of very massive stars from main 
sequence dwarfs to red supergiants will be essayed. Hayashi, 
Hoshi, and Sugimoto (1962) and Stein (1965) have discussed 
the stars of lower mass in similar detail. In particular, 
the evolutionary history of the stellar center, surface, and 
energy-producing regions may be used to summarize a star's 
total evolution. We shall investigate each of these points 

in turn for the case of a star of 30 M . 

O 

jj) Center 

The evolution of the stellar center is shown on a 
log p - log T diagram in Figure 3. It is of particular 
interest to note that the dotted curve of calculated models 
follows closely a line of slope 3, representing Lane's law 
of homologous evolution: 8 c = constant. The major deviations 
occur during the early phases of nuclear burning (both 
hydrogen and helium) when the core expands slightly. The 
slight veering toward isothermal evolution is apparent during 
the gravitational contraction phases. These deviations. 
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however, are small when we compare them with the deviations 
occurring for a star even as massive as 15.6 (Hayashi and 
Cameron 1962, Fig. 6). The reason is the extremely low 
densities and high mass, which minimize effects such as 
electron degeneracy and isothermality . 

b) Surface and Shell Source 

In Paper II an interpretation of the envelope expansion 
was given in terms of the chemical inhomogeneity and the 
U-V plane. We noted that V -* 4 and U-*0asr-*0ina poly- 
tropic envelope of index N = 3. Since the rising luminosity 
in our calculated envelopes centrally condenses the models 
?md increases the mean effective polytropic index toward n = 3 
(or* account of decreasing 8 ) , the U-V curve shifts upward 
and to the left in Figure 4. Thus the increase in total 
stellar radius may be understood from the integration of 
d log r = U - ^ d log q (Hayashi, Hoshi, and Sugimoto 1962): 

Jtocj R = coni ' tanT +■ \ ^ ) + J" "U d J.ocj . 


(36) 
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This relation gives an explanation of the general trend of 

radius expansion beyond the initial main sequence. 

The fluctuations in radius were considered in Paper II 

in terms of the relative importance of radiation pressure. 

1/2 

We showed that since, generally, R ~ (1-8) ' /8T in a poly- 
tropic region, the shell radius and the total radius are 

governed mainly by p and p , respectively, for the core and 

c s 

envelope considered as two separate "polytropes". The results 
for the helium phases are in accord with this picture. 

We shall now discuss further the radius fluctuations 
in terms of the stellar energy requirements and the effect 
of the hydrogen shell source. Therefore we trace the evolu- 
tion from the onset of hydrogen exhaustion to the end of 
helium exhaustion. For completeness, however, we begin with 
a discussion of the initial hydrogen burning phase. 

During this phase of hydrogen depletion in the core, 
the envelope expansion is due to the growing chemical in- 
homogeneity, which reduces U and increases u, (cf. eq. [36]) . 

c 

However, we note that in a completely mixed star, even 
though u monotonically increases with time, the radius 
finally begins to shrink because the decrease in 
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-1.4 

eventually outruns the increase in u c ~ (1+X c ) (cf. e< 3* 

[16] , Paper I) . 

In general, the rise of luminosity during hydrogen 

depletion in a star may be understood by the increase of 

4 

mean molecular weight, through the rough relation L ~ u . 

More precisely, as hydrogen is consumed near the center, 
the energy flux expands the star, and the central tempera- 
ture rises to offset the decreases in density and hydrogencon 
tent. The consequent steepening of the temperature gradient 
must then be maintained by an increasing luminosity, 
according to the radiative transport equation. 

During the hydrogen exhaustion phase, gravitational 

contraction of the core helps to raise T (in accordance 

c 

with the virial theorem) and hence T s . However, since T g and 

the time scale are much too small to displace q by nuclear 

s 

burning, nowhere in the star does |i change significantly. 

Consequently, L remains nearly constant, thereby freezing 

the non-dimensional structure of the envelope. Because 

P. ~ T ^ for such a configuration, the whole radius shrinks, 
c 

This behavior is similar to the evolution during the pre- 


main-sequence contraction phase. 
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During the brief gravitational contraction phase 
following hydrogen exhaustion, the accelerated central 
condensation produces a sudden rise in total luminosity. 

The reason for this is twofold. First, in stars of high 
mass, the large core cannot withstand an energetic gravita- 
tional contraction, so that a temperature gradient persists 
in it. Since this gradient steepens on account of the marked 

increase of T compared with T (because of the sensitivity 
c s 

of L to temperature) , a suitably high luminosity in the core 
H 

is required to support the gradient. Maintenance of such 
luminosity near the center leads to a higher total luminosity 
in the case of a gravitational source of energy, because P 
is more uniform throughout the core than & H . Second, the 
gravitational infall of the shell raises its temperature and 
thereby augments the shell luminosity, too. Hence the whole 
stellar luminosity increases significantly. 

In stars of low mass with isothermal cores less 
than the SchOnberg-Chandrasekhar limit (q = 0.1), the 
luminosity is due almost entirely to shell burning; the 
increase of luminosity will then be due to the temperature 
rise at the shell, caused by slow contraction of the isothermal 
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core. In stars of all masses, however, the radius expansion 
can be ascribed to the marked increase of energy flow into 
the envelope. 

During the early stages of helium burning when the time 

scale is long, nuclear burning in the shell causes a large 

increase of q . To compensate for the outward motion of the 
s 

shell (which tends to reduce T g ) and the depletion of helium 
in the core (which increases 4) , the central temperature 
continues to rise. Although the shell burning is conse- 
quently maintained in order to meet the total flux requirements, 
its luminosity drops gradually because of the slowly decreasing 
temperature. Hence the total luminosity of the star is merely 
kept at a constant value and not increased in accordance with 
L ~ u 4 , as for a point convective source. To keep the shell 
burning at a sufficiently high temperature, the total radius 
must shrink. This may be seen simply from consideration of 
the frozen non-dimensional structure of the envelope (since 
L is nearly constant). As t decreases for larger q^, R must 

also decrease in order to maintain T ~ t R”^. 

s s 

When q g is stabilized in the later stages of helium 
burning (because of the short time scale due to low Y ) , 
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the gradual increase in luminosity expands the envelope. 

This expansion is reflected in the monotonic decrease of L 

H 

(through lower T , which, because of the near constancy 

of t g , implies an increased R) . The deficiency of is met 

by an increase in L caused by the growing mass fraction, 

core 

q^, of the convective core, which mixes down fresh helium. 

This is in contrast with the hydrogen exhaustion phase where 

q„ (and consequently L ) decrease. The total luminosity 
4 core 

of the star now rises more noticeably as q remains nearly 

s 

4 

constant, in accordance with the trend indicated by L ~ u . 
However, since the helium content is already low, jj. and 
therefore L do not change as much as during hydrogen burning. 

During the helium exhaustion phase, the shell source 
supplies the small deficiency in the core luminosity. Thus 
the stellar structure behaves in a way more like the gravi- 
tational contraction phase following hydrogen exhaustion 
rather than the hydrogen exhaustion phase itself, in which 
the shell source and the chemical composition discontinuity 
were just being set up. Analogously, the luminosity rises 
and the radius expands. By the end of helium burning, the 
total luminosity is already so high that any small increase 
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produces a large central condensation of the envelope. This 
explains why, in very massive stars, the luminosity remains 
practically constant throughout helium burning (which resists 
central condensation) , and thereafter produces a tremendous 
envelope expansion with only a moderately small brightening. 
The high luminosity also explains the complicated occurrence 
of convectively unstable regions in the envelope. In such 
convective regions the temperature gradient is practically 
fixed, so that the increasing temperature on the inside 
promotes expansion of the radius. 

£) Convective Core 

As nuclear fuel is ignited in the central regions of 
stars, a convective core is formed if either the temperature 
exponent of the energy generation rate or the radiation 
pressure ( 1— 0) is high enough. This convective core grows 
outward in mass fraction very rapidly at first, reaching a 
maximum extent, whereupon it subsequently retreats. However, 
its detailed behavior is very different if the nuclear energy 
is due to hydrogen or helium burning. In the case of hydrogen 
burning, maximum core extent occurs at the end of the pre- 
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main-sequence contraction, that is, when X is still equal 

c 

to X . However, in the case of helium burning, it occurs at 
e 

the onset of core contraction during helium exhaustion, when 
Y c ^ 0.01. The difference in behavior is essentially caused 
by the electron scattering opacity, and therefore occurs only 
in stars where this opacity dominates. 

The reasons may be seen in detail, as follows. The 
boundary of the convective core is defined as the point 
where (n+1) = (n+1) . From Paper I, we then have that 


it 


o. IS 

16 ire. G-M 


" 1 tX 

_ 1-P 


L (ri + 1 




(37) 


Since the fractional change of (n+1) is small compared 

ad 

with the fractional changes of the other variables in equation 

(37), we may consider (n+1) as constant. Now, during the 

ad 

pre-main-sequence contraction phase, L and 1-0 are approxi- 
mately constant. However, L(r ) rises rapidly due to the 

4 

onset of hydrogen burning near the center, the balance in L 

being supplied by the envelope contraction. Thus q^ ~ L(r^), 

or q increases. 

4 



43 


Since hydrogen burning fully sustains the core during 

the H-phase, the density increase is halted, permitting 1-3 

to rise with temperature. Thus L = L(r ) ~ 1-8 (roughly), 

4 4 

with the result that q, ~ h ~ 1 + X . Hence during hydrogen 

444 

depletion, the opacity decreases and the core mass shrinks. 
These considerations supplement those of Sakashita, Ono, and 
Hayashi (1959) , who showed that the convective core could at 
least not grow in mass. 

Helium ignition in the core behaves in a similar way to 

hydrogen ignition, resulting in an increase of q . However, 

4 

during the He-phase, the now constant opacity (*. = 0.19) does 

not regulate the core size as it did during the H-phase. 

Since the structure is altered by the presence of composition 

discontinuities and a shell source, it is also no longer true 

that L(r ) ~ 1-8 . In fact, L(r ) increases by a much larger 
4 4 4 

factor than 1-8 can. Hence q ~ L(r )/(l-B ) , or q increases 

4 4 4 4 4 

as helium is depleted. 

During helium exhaustion, gravitational contraction in 

the core begins to replace the nuclear source of energy. 

Since the gravitational source is uniform and distributed 

well beyond the convective region, L(r ) decreases. The 

4 

density rises considerably and checks the increase of 1-8 , 

4 
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which, as during all gravitational contraction phases, 
remains roughly constant. Thus it is only at the end of 
helium burning that the convective core begins to shrink, 
on account of the decreasing luminosity at q^. 

XIII. LATE EVOLUTIONARY PHASES 

Whether the high-temperature neutrino processes exist 
in nature or not, their inclusion in our models up to the 
end of helium exhaustion does not substantially jeopardize 
the results since the neutrino luminosity rises to only a 
few per cent of the photon luminosity at this stage. To 
calculate further a rough time scale of the ensuing evolu- 
tionary phases, we must discuss separately the cases of 
evolution with and without neutrinos. In both cases we 
assume that the initial core composition is pure oxygen. 

a) Photonic Evolution 

Under normal evolutionary conditions, gravitational 
contraction will heat up the core until oxygen burning 
sustains the structure. We assume that the gravitational 
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energy release, which is taken to be uniform throughout 
the core, may be represented by 

. 3 y]L _ U_ 

e 3 = 2 3r iMfr) M ' ( 38 > 

Solving for the time, we have 


3 A. M ^ 1 e. o 

T L * 


(39) 


Equation (39) will no longer be valid when L > L . 

nuc g 

This will happen when t = t . If we assume a nuclear 

nuc g 

energy generation rate in the form of equation (17) , then 
it may be shown for a polytrope of index 3 that (Reeves 
1963) 


£/e c = (u +6:o * r ) (40) 

g 

With q^ - 0.4, w^ = 0.5, L = 10 L^, and the evolutionary 

3 9 r> 

relation p ^ T , we obtain T = 1.5 x 10 °K when 

c 

f^nuc = C . By using this final temperature and an initial 

8 

temperature of 3.7 x 10 °K, equation (39) yields 

4 

t — 3 x 10 years. 
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Thereafter, oxygen burning and successive gravita- 
tional contractions interrupted by still more advanced 
nuclear burning phases up to iron formation supply the 
stellar luminosity. Since the gravitational phases will 
be negligibly short compared with the nuclear phases, 
we consider only 


r- EZM^/Lu.^. (4i) 

Oxygen burning will be the dominant nuclear phase. We 

estimate roughly that Z = 1, = 0.3, w^ = 0.5, and 

0 1 7 

L = 10 L^. With E = 9 x 10 erg/gm for the conversion 

of oxygen into iron, we obtain t = 3 x 10^ years. Even 
if the initial core composition were pure carbon, the 
lifetime of the star as a red supergiant would not be 
substantially increased. 

b) Neutrinic Evolution 

On the assumption of neutrino emission, the longest 
of the late evolutionary phases will be the gravitational 
contraction following helium exhaustion, when the evolution 
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is at first primarily limited by photon radiation. Thus 

equation (39) will apply up to a certain temperature at 

which 6 = G . This temperature may be calculated as 

v 9 

above. With w =1 and L = 4 x 10 5 L , we obtain 
g © 

T = 6.1 x 108 °K. Hence t = 7 x 10 3 years, 
c 

g 

At a temperature around 6 x 10 °K, three changes in 

the stellar luminosity will take place. First, neutrino 
emission will predominate over the photon radiation. 

Second, photoneutrinos will be replaced by pair annihila- 
tion neutrinos as the dominant neutrino source. Third, 
the helium shell will be ignited. However, since the total 
photon luminosity will not rise very much in comparison 
with the increasingly powerful neutrino luminosity (Reeves 
1963), the helium shell may be neglected. 

The evolution of the star is thus limited only by the 
rate of neutrino emission from pair annihilation, that is, 

3 

= L g* Under the assumption that p « T , the time scale 
is given by (Stothers 1963a) 



This yields t = 6 x 10 3 years. 


( 42 ) 
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Oxygen burning may briefly sustain the neutrino 

9 

losses at a temperature around 2 x 10 °K. For a poly- 

tropic core mass of 9 M^, Reeves (1963) estimates the 
oxygen burning lifetime at 2 months or less. Subsequently, 
gravitational contraction must supply the neutrino losses 
while further nucleosynthesis up to iron formation takes 
place in the core, on a time scale of minutes. Thus the 

total lifetime from the end of helium exhaustion will be 
4 

about 1.3 x 10 years, on the neutrino hypothesis. 


XIV. PRE- MAIN- SEQUENCE CONTRACTION 

Chandrasekhar (1939) has given an expression for the 
pre-main-sequence contraction time: 


3y-M- GJA 2 

Y ~ 1 lr' * 


(43) 


If we assume for the ratio of specific heats y = 5/3 and 
for the polytropic index N = 3, then equation (43) yields 

4 

2.3 x 10 years. The implicit assumption that L = constant 
is a good approximation for very massive stars. In fact, t 
as deduced from equation (43) is reasonably close to values 
for 30 Mq interpolated from the approximate models of Tanaka 
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and Sakashita (1965) and from the more exact models of Ezer 
and Cameron (1965) . 

XV. FINAL REMARKS 

The time scale of evolution for the different phases 
of a 30 star is given in Table 5. The helium burning 
lifetime is one-tenth of the hydrogen burning lifetime? this 
fraction is somewhat larger than for 15.6 (Hayashi and 
Cameron 1962) because of the higher initial hydrogen abun- 
dance assumed there. The other phases contribute negligibly 
to the total lifetime of slightly less than 6 million years, 
except in the case of no neutrino emission. In that case, 
oxygen burning to iron lasts a period comparable with the 
helium burning lifetime. 

The evolution of the star has been plotted on the theore 
tical H-R diagram in Figure 5. The shape of the track is 
very similar to the one computed for 15.6 , except that 

at 30 the luminosity during helium burning increases very 
little . 

The star sweeps out nearly the whole B spectral range 
during helium burning as a supergiant, from B8 to Bl. By 
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the end of helium exhaustion, it passes quickly across the 
Hertzsprung Gap and becomes a K supergiant. Further gravi 
tational contraction in the core is expected to bring the 
star into the region of the brightest M supergiants. 
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Distribution of Chemical Composition in a Star of 
30 M during Helium Burning 
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TABLE 2 


Evolutionary Models of a Star of 30 M 0 during Helium Ignition (I) Phase* 


Models 


b 

1 

2 i 

3 ! 

i 

4 i 

5 

6 

f 

log C 

-3.145 

1 

-3.140 j 

\ 

-3.137 | 

-3.135 ? 

-3.133 

-3.131 

*1 

0.772 

0.774 

0.775 : 

0.775 

0.776 

0.777 

q 2 q 3 } 

0.461 

0.456 ; 

0.453 - 

0.451 j 

0.449 

0.447 

B s 1 

0.455 

0.448 j 

0.444 j 

j 

0.441 j 
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0.436 

lo 9 T s | 
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i 

1 

7.722 j 

s 
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t 
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0.669 I 

i 

0.674 | 

J 

1 
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0.687 

log (r g /R) 

-1.968 

i 

-2.099 ! 

1 

1 

-2.185 

-2.247 1 

j 

-2.312 ■ 

! 

-2.366 

Us : 

0.288 

0.260 1 

0.245 

0.235 i 

i 

0.226 

J 

j 

0.219 

i 

v s ! 

3.943 

3 . 984 

4.008 

1 

4.022 | 

} 

1 

4.035 
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(n + l) s 

3.978 

3.981 

3.981 

3.983 j 
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3.984 

q 4 

0.033 
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0.150 | 
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0.210 

t 

t 

B c ; 

0.695 
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0.676 
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0.672 

i 

log T i 
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8.243 

8.256 i 
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8.265 ' 
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8.269 

log o c 
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2.716 i 
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J 
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2.733 

V L 
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0.535 j 

0 . 566 
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i 0.629 
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W* 
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i 

0.147 ! 

1 
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| 0.297 
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L g /L 

0.442 
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0.224 
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| 0.074 


log (L/L 0 ) 

5.553 

i 

5.558 ! 
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5.563 

j 5.565 

1 ; 
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log (R/R 0 ) 
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2.139 

log T e 
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4.163 

4.137 
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4.086 

4 

t (10 years) 

0.00 

0.24 

0.42 

0.58 

0.84 

1.50 


*Also, q = 0.328; X = 0.018; \ = 1.552; and Y = 0.97, X - 

s s s ^ 

convective core. 
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Evolutionary Models of a Star of 30 M during Helium Depletion (He) Phase 
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TABLE 4 


Evolutionary Models of a Star of 30 M Q during Helium 
Exhaustion (X) Phase* 



Models 
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2 

3 

4 
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\ 
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FIGURE LEGENDS 


Fig. 1. — Evolution of the structural zones of a star of 
30 M q during the phases of helium ignition, 
depletion, and exhaustion. 

Fig. 2a,b. — Distribution of hydrogen and helium as a 

function of mass fraction during the helium deple- 
tion (He) and exhaustion (X) phases. Dots represent 
zonal boundaries. 

Fig. 3. — Evolution of the center of a star of 30 M from 

© 

the initial main sequence into helium exhaustion 
(Models H0-X2). Circles represent the first stages 
of the various phases. The line has a slope of 3. 

Fig. 4. — Evolution of a star of 30 M^ in the U-V plane, 

during the helium depletion (He) and exhaustion (X) 
phases. Dots represent zonal boundaries. Dashed 
lines refer to discontinuities in the chemical com- 
position. Model X5 is based on the assumption that 



Fig. 5. — Evolution of a star of 30 M on a plot of lumi- 

© 

nosity versus effective temperature, from the initial 
main sequence through helium exhaustion. Circles 
represent the first stage of the various phases. 



